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Abstract 



We show that the mathematical structure of Gibbsian thermody- 
namics flows from the following simple elements: the state space of a 
thermodynamical substance is a measure space Q together with two 
orderings (corresponding to "warmer than" and "adiabatically accessi- 
ble from" ) which satisfy certain plausible physical axioms and an area 
condition which was introduced by Paul Samuelson. We show how 
the basic identities of thermodynamics, in particular the Maxwell re- 
lations, follow and so the existence of energy, free energy, enthalpy and 
the Gibbs potential functions. We also discuss some questions which 
we have not found dealt with in the literature, such as the amount of 
information required to reconstruct the equations of state of a sub- 
stance and a systematic approach to thermodynamical identities. We 
illustrate the power of these methods by giving in detail explicit com- 
putations for several real gases which, in the textbooks, are normally 
only obtained for the simple case of an ideal gas. 

This article is dedicated to the memory of Paul A. Samuelson, who had 
the initial insight which led to our approach and who accompanied its devel- 
opment with a barrage of new questions, ideas and encouragement. 

1 Introduction 

The subject of thermodynamics is notoriously difficult for mathematicians. 
V.I. Arnold [Ar] famously put it in a nutshell as follows: 

Every mathematician knows that it is impossible to understand 
any elementary course in thermodynamics. 

He continues by explaining that 

the reason is that [the] thermodynamics is based on a rather 
complicated mathematical theory, on [the] contact geometry. 

It is the purpose of this note to present an axiomatisation which is mathe- 
matically transparent, avoids anthropomorphisms, is elementary (no contact 
geometry) and preserves all of the structure of the classical theory. It also 
allows us to carry out exphcit computations for the classical cases (ideal gas, 
van der Waals gas) in a simple and general way which works for virtually any 
of the standard models for real gases. We illustrate this for an equation of 
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state suggested by Feynman to allow for the fact that, for real gases, the adi- 
abatic index is far from being constant and for a new model which combines 
the advantages of the van der Waals and the Feynman gas. 

Wc also discuss a theme which seems to us to be of eminent practical 
significance but which we have never seen treated in the literature — namely, 
how much information is required to recreate the equations of state of a 
substance. Thus we show that knowledge of all of the isotherms and just two 
adiabats (dually two isotherms and the adiabats) suffices. 

Our system of axioms has two basic ingredients — firstly a measure space, 
i.e., a set with a a-finite positive measure (which describes the possible states 
of a thermodynamic system), provided with two pre-orderings, "warmer 
than" and "adiabatically accessible from". We show that if these order- 
ings satisfy several physically plausible conditions, then they are induced 
by numerical functions — empirical temperature and empirical entropy in the 
thermodynamical case. 

We then show that an area condition which was made explicit by the 
distinguished economist P. A. Samuelson ensures the existence of essentially 
unique absolute temperature and entropy. These have the property that the 
four Maxwell relations then hold. If we then interpret these as integrability 
conditions, we can deduce the existence of the four energy type functions 
of thermodynamics. This is in contrast to many standard treatments where 
these relations are deduced from the existence of the energy functions via 
the Schwarz Lemma. In a final section, we present a unified and systematic 
approach to the classical thermodynamical identities between derived quan- 
tities and discuss some known and new models of real gases in much more 
detail than one finds in standard treatments. 

2 Samuelson's vision 

Samuelson noted that classical thermodynamics and economics are related 
by a common search for an optimising basis for observed behaviour. In 
thermodynamics the observed isotherms and adiabats are hypothesised to be 
derived from the minimisation of a scalar quantity "energy" . In economics, 
the observed input demand functions are hypothesised to be derived from 
the maximisation of the scalar quantity "profits" . Deriving a test for these 
hypotheses then becomes a common task of both disciplines. 

Thus we can interpret results of Maxwell as establishing the equivalence of 
the existence of an energy function which is minimised with the fact that the 
isotherms and the adiabatics fulfill a simple and natural geometric condition 
which we called in an earlier paper the "S'-condition" — see below for a precise 
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formulation of the latter. 

Such a condition is certainly implicit in Maxwell's argument. However, it 
is not stated explicitly. Samuelson claims no priority for having noticed via 
this diagram that constrained optimisation implies this relationship but we 
have been unable to find another reference stating the equilibrium condition 
in such a geometrically simple way. Thus although the above area condition 
appears implicitly in many areas, we have found no previous instance of its 
explicit formulation in the literature. 

In [Col] we obtained a number of equivalent formulations of this condi- 
tion, notably a rather lengthy partial differential equation in two arbitrarily 
chosen functions which have the given curves as contours (corresponding to 
empirical temperature and entropy in the thermodynamical situation). We 
then showed that when this equation is satisfied there are canonical recal- 
ibrations of empirical temperature and entropy for which the Jacobian is 
identically "1" and this implies the existence of an energy function and the 
vahdity of the Maxwell relations for the recalibrated quantities. Note that, as 
in Maxwell's original treatment, we derive his relations from an area condi- 
tion in {p, y)-space, the existence of recalibrations with the J = 1 condition 
being a consequence of the S condition. This paper can thus be viewed as 
an attempt to rehabilitate and perhaps clarify the Maxwell approach. 

We shall use this theory to discuss existence and uniqueness of solutions 
to the above partial differential equation. In particular, we show that given 
one family of curves (say the isotherms) and two other members of the other 
family, then this uniquely determines the family of adiabats (for a more 
precise mathematical formulation, see below). 

The Maxwell/Samuelson area condition thus establishes an important 
duality between isotherms and adiabats and, more generally, between any two 
suitable families of level curves which are derived by minimizing energy in two 
different constraint regimes. These facts can be used to show how to exploit 
this duality to derive some explicit formulae for the dual functions (adiabats 
given isotherms, isotherms given adiabats) for two standard textbook cases, 
the ideal gas and the van der Waals gas. 

Recently there has been renewed interest in the derivation of the Maxwell 

dip, V) 

relations in thermodynamics from the Jacobian identity = 1 (see, for 

o[T, S) 

example, [Ri]). As is well known, this identity means that the the correspond- 
ing map from the (T, S')-plane into the (p, y)-plane is area preserving, and so 
this approach links thermodynamics to such geometrical mechanics 

where area-preserving mappings play a central role. Maxwell was well aware 
that his relations had a geometric foundation and indeed he used Euclidean 
geometry "to get his four identities in an amazingly obscure way" [Am] . Per- 
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haps because of its opacity, Maxwell's use of geometry to derive his relation- 
ships seems to have largely disappeared from the physics textbooks. In the 
discipline of economics, however, exactly this obscure geometrical argument 
of Maxwell was the subject of the acceptance speech of Paul Samuelson, the 
first American Nobel laureate in Economics. The Maxwell/ Samuelson area 
condition thus establishes an important duality between isotherms and adia- 
bats and, more generally, between any two families of level curves which are 
derived by minimizing energy in two different constraint regimes. 

3 Orderings and the axiomatics 

In this section we discuss briefly a topic which is probably the basic problem 
of the theory of measurement — when is a physical quantity described (in a 
meaningful fashion) by a number? If this can be done, then the values of the 
quantity in question can be compared (as in "warmer than", "worth more 
than" etc.) and the basic question is when the converse holds , i.e., when is 
such an ordering induced by a numerical function ("temperature", "price")? 
This has been examined and re-examined countless times (to our knowledge, 
the first rigorous formulation of a mathematical theorem of this sort was 
due to Dcbreu [De] who gave a sufficient condition for a preordering to be 
induced by a utility function). We give a brief discussion, firstly for the sake 
of completeness and secondly because we would like to emphasise what we 
regard as the central point, namely that the real fine has a simple and elegant 
characterisation as an ordered space. The latter fact is wellknown (although 
perhaps not quite as wellknown as it should be) but we have never seen it 
used on the problem we are now addressing. The opening pages of Maxwell's 
treatise [Ma] give a lucid treatment of temperature as an ordering. 

3.1 Orderings and utility functions 

Suppose that we are given a set Q and a surjective mapping / from it onto 
the real line. (Since we are only concerned with the order-theoretic aspects of 
the latter, we can, at will, replace it by any order-isomorphic set, for example 
an open interval, in particular the half-line). Then / induces a preordering 
<f where we define: x <f y ii and only if f{x) < f{y)- This preordering has 
the following properties: 

1. it is total and has no largest or smallest element; 

2. it is order complete; 
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3. Q contains a countable, order dense subset. 

For our purposes it will be convenient to restate these properties in terms of 
the family of predecessor sets of elements. Thus, we define, for each a e Q 
the sets 

Aa — {x G fl : X <f a}, Ua = {x E il, : x <f a}, 
which have the following properties: 

1. The Aa are distinct and totally ordered by inclusion; 

2. the countable subfamily {Aq : q e Q} is order dense; 

3. U^a = ^^, n^a = 0; 

4. The family of the A^ is closed under intersections and for each a, 

Aa = f]{Ap ■./3>a}. 

The family of the C/'s satisfies the corresponding properties, except that 4. 
is replaced by 

4'. the family is closed under unions and for each a, Ua — U{^/3 ■ P < o;}- 

Further, if a < /5 < 7, then 

Aa (Z U/3 (Z A/3 (Z (Z A^. 

(Note that if we start with a family Aa as above, and define Ua — {JiAp : 
P < a}, then the above conditions are fulfilled). 

We remark at this point that in this paper the inclusion C will always be 
exclusive, i.e., A G B implies that A and B are distinct. 
We are interested in the following converse statements: 

Theorem 1 Suppose that we are given a family A of subsets of Q which 
is totally ordered by inclusion, is closed under arbitrary intersections and 
satisfies the properties 

1. if Ae A, then A = [\{B eA,Ac. B}; 

2. there is a countable subset Aq which is order dense in A; 

3. {^A = % and{jA = n. 
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Then there is a surjective mapping f from Q onto R such that A is the family 

{f < a: a e R}. 

The simple proof of this result follows from the following standard order- 
theoretical characterisation of the real line: 

Theorem 2 Suppose that we have a set A with a total ordering such that 

1. A has neither a smallest nor a greatest element; 

2. A has a countable, order- dense subset; 

3. A is order complete. 

Then A is order-theoretically isomorphic to the real numbers. 

This in turn follows easily from the following characterisation of the rationals, 
which is due to Cantor: 

Theorem 3 Suppose that we have a set A with a total ordering so that 

1. A has neither a smallest nor a greatest element; 

2. A is countable; 

3. A is dense in itself (i.e., if x < z in A, then there is a y E A with 
X < y < z). 

Then A is order-theoretically isomorphic to the rational numbers. 

The idea behind the proof of theorem 1 is now simple. We introduce the 
following equivalence relationship on Jl: a; ~ y if and only if x and y are in 
exactly the same sets oi A e A, i.e., for each A, x E A ii and only ii y E A. 
Then the quotient space fl\^ has a natural order structure which satisfies the 
properties which characterise the real line. The required mapping / is then 
the canonical one from fl onto the quotient space. 

An important point is that the / in theorem 1 is not uniquely determined. 
We can replace it by any F — 0o/ where (p is an arbitrary order isomorphism 
of the line. We call such an F a recalibration of /. In the general situation 
there is no canonical choice of F. This will be crucial in the following. 

There are two refinements of theorem 1 which will be of particular interest 
to us. Firstly, if fl is provided with a suitable a algebra and each A e A is 
measurable (i.e., a member of the algebra — we then say that the ordering is 
measurable) , f will be measurable. Secondly, if is a topological space and 
each A E Ais closed and, further, for each A E A, U = [J{B E A : B G A} 
is open, then / is continuous. There are corresponding conditions which 
ensure semi-continuity. 
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3.2 The axiomatics 



We are now in a position to state the four axioms which describe the math- 
ematical structure of a thermodynamical theory: 

1. The states of a thermodynamical system are specified by the points of 
a set with a positive cr-finite measure /i. 

2. fl is provided with two families .Atcmp and ^cnt of measurable subsets 
which satisfy the conditions of theorem 1 (and so the preorderings are 
induced by numerical functions which we denote by t and s). 

3. for each A_i (Z Aq G Ai in ^temp and B_i (Z Bq (Z Bi in ^ent we have 

\ ^o) n (B, \ Bo))iJi{Ao \ ^-i) n {Bo \ 
= //((Ai \ Ao) n {Bo \ B_,))^{{Ao \ n (Si \ Bo)). 

4. Condition 3. means that if t and s are the (measurable) functions 
which induce these orderings, then the image measure of /x in under 
the mapping u h-> (t(a;),s(w)) splits multiplicatively. We further as- 
sume that this measure is equivalent to Lebesgue measure on the plane 
(equivalent in the sense of being mutually absolutely continuous) . 

It follows from this that for each A G .Atcmp) and for each A G .Acnti 
li{A \ \J{B e A: B C A}) = 0, and that for each pair A G Ai from A and 
B G Bi from B wc have 

f,{{A^\A)n{B^\B)) 

is strictly positive and finite. 

We remark that these axioms are physically natural and plausible. (In 
classical thermodynamics, Q is {p, y)-spacc and the measure is interpreted as 
mechanical work). Since this article is, despite its title, one in mathematics 
rather than in physics, we will not go into this in detail; but we emphasise 
that the area condition, in particular, is not a deus ex machina inserted to 
save the day but has a natural physical justification. A further point, which is 
important more for philosophical reasons, is that they do not explicitly refer 
to the real numbers (compare the axioms for Euclidian geometry, particularly 
in the form as perfected by Hilbert). 
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3.3 The existence of absolute temperature and entropy 



If we only assume 1. of the above axioms, s and t are, as is the function / in 
the general result, not uniquely determined. A crucial point of our treatment 
is that in the presence of conditions 3. and 4. above, then there are (essen- 
tially) unique such choices, which we call the canonical recalibrations. For 
the condition 4. means that the image measure has the form a{u)b{v) dudv 
where we denote the coordinates in by (u, v) and a and h are locally 
Lebesgue-integrable functions whose reciprocals are also locally integrable. 
Now if we replace the two functions t and s by the recalibrations T — (f) ot 

and S — ip o s, where is a primitive of - and of -, then we obtain the 

a b 

following result: 

Theorem 4 Suppose that the above axioms are satisfied. Then we can choose 
the functions T and S so that the mapping cu i— )■ {T{u), S{u!)) is area- 
preserving. 

This choice is unique up to suitable affine transformations (loosely speak- 
ing, we can choose the zero point and a change of scale — c.f. the difference 
between the Celsius and Fahrenheit systems). 

In thermodynamics, these canonical calibrations arc called absolute tem- 
perature and entropy (as opposed to empirical temperature and entropy). 

We now proceed to show that these axioms imply the usual contents of 
elementary treatments of thermodynamics. 

As we shall see shortly, this choice of calibration is crucial since the fact 
that the area condition holds is equivalent to each (and hence all) of the 
four Maxwell relations. Since the latter can be interpreted as intcgrability 
conditions, they ensure the existence of the four energy type functions of 
thermodynamics. (Once again, these are purely mathematical facts, but 
the underlying motivation from thermodynamics is the principle of Joule- 
Maxwell on the mechanical equivalence of heat) . 

4 Samuelson configurations 

4.1 The area condition 

We emphasise at this point that up till now smoothness (except in the very 
mild form of measurability) has played no part in our considerations, neither 
in the formulation of the axioms nor in the derivation of the canonical recal- 
ibrations. This is as it should be, for philosophical reasons but also because 
the presence of phase transitions makes it clear that in real substances we 
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can and should expect the isotherms and adiabats to have "corners". We 
now turn to the case where they are smooth, i.e., we suppose that fl is the 
plane or some suitable subspace (in thermodynamics usually the posi- 
tive quadrant) and that the functions which induce the ordering are smooth 
(in the sense of being infinitely differentiable) as are their level curves. It is 
convenient to use the mathematically neutral notation x and y for the coordi- 
nates in the plane which we shall initially regard as the independent variables 
and u and v for the two (potential) functions. As shown in [Col], the area 
condition is equivalent to the fact that the functions u and v satisfy a certain 
non-linear partial differential equation of the third order which is displayed 
exphcitly there. Familes of level curves satisfying this partial differential 
equation are thus of great interest in the study of optimizing systems. 

In the following, we concentrate on the two foliations consisting of the 
level curves of u and v, i.e., the isotherms and adiabatics in the thermody- 
namical context. The area condition can then be formulated as follows. The 
plane (or a suitable part thereof) is foliated by two families of curves — the 
level curves of two potential functions u and v. We assume that these are 
transversal at each point, i.e., the Jacobian J = {uxVy — UyVx) never vanishes. 
Since the regions we consider arc connected (in the topological sense), J can- 
not change sign. Hence there is no essential loss of generality if we assume 
that it is always strictly positive. 

Locally the families of curves form a network which is topologically equiv- 
alent to the standard network of the plane induced by the parallels to the x 
and y axes (i.e., the case where u{x, y) = x and v{x, y) = y). 

We say that the foliations satisfy condition S (or that v is S'-transversal 
to u or that the v-curves are S'-transversal to the u curves) if the following 
holds: for any choice of values c_i < Cq < Ci, and d^i < do < di respectively, 
we have 

area A/area B = area C/area D 

where 

A = {{x,y) : c^i <u{x,y) < Co,do <v{x,y) < di}, 
B = {{x,y) : Cq <u{x,y) < ci.do < v{x,y) < di}, 
C = {{x,y) : c^i <u{x,y) < Co,d_i < v{x,y) < do}, 
D = {{x, y):co< u{x, y) < Ci, d^i < v{x, y) < do}. 

In order to avoid topological problems, we assume that the values of the 
c's and d's are sufficiently close for the above condition on the network to 
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be satisfied. This means that the condition we are considering is a local 
one, as it should be if it is to be equivalent to a partial differential equation. 
However, it is easy to obtain a global form from the local one. 

We refer to the configuration consisting of two foliations which are S- 
transversal as a Samuelson configuration. The S'-configuration consisting of 
the adiabats and isotherms of the ideal gas is one of the most iconic images 
of modern science. 

The precise relation of our result to this question will be made more 
explicit below. For obvious reasons, we discuss the case of the isotherms for 
an ideal gas and for a van der Waals gas in some detail. In particular, we 
show that if a function v is S'-transversal to the function u — xy (i.e., the 
potential defining the isotherms of an ideal gas) and two of v's level curves 
have the form xy'^ = constant for the same 7, then all of them have this 
form, i.e., the adiabatics are precisely those for the ideal gas with exponent 
7- 

4.2 Thermodynamical notation 

Although this article is one on mathematics and not on physics, its main 
motivation comes, of course, from thermodynamics. For this reason, we re- 
call the standard notation and concepts from classical thermodynamics for 
the reader's convenience. The coordinates x and y in the neutral notation 
correspond to p and V in Gibbsian thermodynamics. The choice of the 
latter as independent variables is natural since these two quantities can be 
directly measured. Also the natural meaure on this space (mathematically 
speaking, two-dimensional Lebesgue measure) has a natural physical inter- 
pretation (mechanical work). 

Our starting point is the situation where we are given the temperature T 
and the entropy s as functions of the pressure p and the volume V. We use 
the lower case s to indicate that this is empirical entropy. Absolute entropy 
will be denoted by S. (The standard models do not require a recahbration of 
temperature so that there will be no need at this point to distinguish beween 
lower and upper case T — however, we consider below an interesting model 
due to Feynman where we shall require such a recahbration). We use the 
following dictionary to jump between the purely mathematical notation and 
the thermodynamical one: u corresponds to T, v to s, p to x and V to y. 
For example, the thermodynamical equations 

T = pV, s = pV^ 
of the ideal gas corresponds to 
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u{x,y)^xy, v{x,y)^xy'^. 
For reasons which will be clear shortly, we use the recalibrated form 

u{x,y)^xy, v{x,y) ^ ^—-{lnx + -f\ny) 

7-1 

of these equations. 

It is a consequence of the Maxwell relations that one can define the follow- 
ing four energy type functions (whose definitions we repeat for the readers' 
orientation. They can be found in any textbook on thermodynamics, e.g. 
[La]). Firstly, the energy E = E{S, V) is a function of entropy and volume. 
From this one derives the quantities T (temperature) and p (pressure) by the 
equations 

^E{s^n dE{s,v) 

dS ' ^ dV ' 

Analogously, one has the enthalpy H — H[p,S), from which one derives 
the quantities 

rj. _ dHjp, S) dHjp, S) 

dS ' dp ' 

the free energy F — F{T, V), from which one gets 

dF{T,V) dF{T,V)) 
dT ' ^ dV ' 

and the free enthalpy G = G{p,T), which gives 

dGjp, T) dG{p^ 
^~ &r ' ^ ~ dp ■ 

In the German- language literature, e.g. [La], one employs $ for G and W 
for H. 

We emphasise that in our treatment the logical development is reversed — 
the existence of such functions is a consequence of our axiom system, since 
it follows from the Maxwell relations which in turn are equivalent to the 
vahdity of the ySi-condition (see below) . 

4.3 Canonical recalibrations 

We saw above that if the level curves of u and v satisfy the 5* condition, then 
we can find (essentially unique) recalibrations U — (f)ou and V — '4'ov (where 
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(j) and -0 are diffeomorphisms between, say, intervals of the real line), so that 
the Jacobian is identically one. If we assume that u and v are so calibrated, 
then this means that the diffeomorphism i— )■ {u{x,y),v{x,y)) of the 
plane (or a suitable subset thereof) is area preserving. In this case we say that 
the functions u and v satisfy the S'l-condition, or that v is S'l-transversal 
to u. (The recahbration for the ideal gas which was used above arose in this 
way). 

We now come to the crucial point in our argument. If we write the basic 
equations u = f{x,y), v = g{x,y) in differential form, i.e., as 

du — fidx + /2 dy, dv — gidx + g2 dy 

ifi, h ^'^^ the partials with respect to x, y etc.), we can solve for du and dx, 
say, to get 

du — — dv dy, dx — — dv — — dy, 

9i 9i 9i 9i 

where J is the Jacobi-determinant fig2 — f29i, and so we see that the condi- 

Ou Ox 

tion J = 1 is equivalent to the Maxwell relation — = — — which is an 

oy V ov y 

integrability condition and ensures the existence of a function h of the two 
variables y and v such that u and v are the solutions of the equations 

X- f{y,v)^Q, u- g{y,v) ^0, 

dh dh 
where f(y, v) — — — and g(y, v) = —. (We are using the standard conven- 
oy ov 

du 

tions employed in thermodynamics — thus — denotes the partial derivative 

dy V 

of u, regarded as a function of y and v, with respect to y). 

The proof of this result employs the inverse function theorem and so the 
precise statement is local. The same remark applies to many of the following 
enunciations. 

We shall call such a function h a geometric energy function since in certain 
situations where the foliations arise as the level curves of suitable physical 
quantities it corresponds to the energy of a system. However, in such situa- 
tions, the energy function satisfies some structural properties (monotonicity, 
convexity) which have natural physical interpretations, and these are of no 
direct relevance in our considerations below. In a similar manner, we will 
talk of geometrical adiabatics associated with famihes of isotherms, or geo- 
metrical isotherms associated with families of adiabatics respectively. 

The existence of the above energy functions is one of those facts which 
have been discovered and rediscovered time and again in the history of math- 
ematics. We have traced it as far back as to Gaufi [Ga] who used it to describe 
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all equivalent projections (in the sense of mathematical cartography) and it 
appears in contact geometry (under the name of a generating function). Of 
course, as remarked above, it has long been used in thermodynamics. 

It follows from the above observation that we have a remarkable symmetry 
(corresponding to the Maxwell relations in thermodynamics). If we start with 
a given energy function, we can as above calculate u and v as functions of 
X and y. The energy function arises from the process of replacing x and y 
as independent variables by y and v. There are four such possibilities (the 
interesting ones for us are those with y and v, x and v, y and u and x and 
u respectively as independent variables), each of which is associated with an 
"energy function" (as we noted above, in thermodynamics they are called 
energy, free energy, enthalpy and free enthalpy respectively). Hence any 
one such function automatically defines three others. (In fact, the situation 
is more complicated than described here. This is due to the fact that we 
are relying on global solvability of the corresponding non-linear equations. 
The general results we use employ the inverse function theorem and so only 
guarantee local solubility. In many concrete situations which we compute, 
wc do, of course, have global invertibility and hence the kind of symmetry 
evoked here). 

In the case of an ideal gas, the permutations of the various variables can 
be computed by hand and are valid globally — we include the formulae below 
for completeness. We have also added a more general case since it displays 
the fact that the familiar presence of the logarithm in the expression for the 
entropy of an ideal gas is in a certain sense unique to this case. Already 
the van der Waals gas offers difficulties here and we shall shortly develop an 
alternative method of computing these energy functions which is often more 
practical and doesn't require us to compute these permutations. 

5 Thermo dynamical identities — an anthologie 

misonnee 

5.1 The basic machinery 

We suppose that u and v are given as functions / and g oix and y. Thus u — 
f{x,y), V — g{x,y) and, when J — 1, simple manipulations with differential 
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forms proved the basic identities: 
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-9i 


du 


+ 


Jidv 


du — 


^dv 

92 


+ 


j_ 

92 




dy = 


j_ 

92 


dv 




s^dx 

92 


du — 


^dv 

91 




±_ 

91 


dy, 


dx — 


J_ 

91 


dv 




a^dy 

91 


dv — 


fdu 

Jl 


+ 


1 

fl 


dy, 


dx — 


1 

fl 


du 




^dy 


dv — 


fdu 

h 




1 

h 


dx, 


dy = 


1 

/2 


du 




T" dx 



where we have highhghted the expressions which correspond to the Maxwell 
relations. 

In thermodynamic notation these are 



dT = fidp + f2dV, dS 
dp = 92 dT — g2 dS, dV 

dT = ^dS+—dp, dV 

92 92 

dT - ^^5- dV, dp 

9i 9i 

dS = ^dT+^dV, dp 
Jl Jl 

dS = ^dT-^dp, dV 

J2 J2 



91 dp + 92 dV 
-gi dT + f\dV 

— dS — — dp 

92 92 

— dS-^dV 
91 91 

— dT-^—dV 

h 92 

dT-fdp 

J2 J2 



where we are using the key: u^T,v-^S,x-^p, y-^V introduced above. 
In order to isolate the underlying patterns, we now use a numerical code. 
Thus 

u 3^ T 
v ^ S 

X — > 1 ^ p 
y ^ 2^ V. 

Partial derivatives will be denoted by triples in brackets. (3,1,2), for exam- 
pie, denotes fl, in the neutral notation, in the thermodynamical one. 
In general, k) denotes the partial derivative of variable i, regarded as a 
function of the j-th and fc-th variable, with respect to the j-th variable. 
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By reading off from the above list, we can express each partial derivative 
of the form k) in terms of fi, f2, Qi, 92 as follows: 



(3,l,2) = /i, (3,2,l) = /2, (4, 1,2) = 51, (4,2,1) =52; 



[1,3,4) = g2, (2,3,4) = -51, (1,4,3) = -/2, (2,4,3) = A; 



(3,4,1) = ^, (2,4,1) = 1, (3,1,4) = 1, (2,1,4) = -^; 

92 92 92 92 



XI f 
(4,3,l) = f, (2,3,1) = -, (4,1,3) = -, (2,1,3) = -^; 

J2 J2 J2 J2 



(3,4,2) = :^, (1,4,2) = i, (3,2,4) = -i, (1,2,4) = -^; 

9i 9i 9i 9i 



(4,3,2) = ^, (1,3,2) = |, (4,2,3) = |, (1,2,3) = -^. 

Ji Ji Ji Ji 

We can then express any derivative (a, b, c) in terms of ones of the form 
{d, 1,2) or (e,2, 1). Thus the four derivatives with x and v as independent 
variables are as follows: 



(3,4,1) 
(2,4,1) 
(3,1,4) 

(2,1,4) 



(3,2,1) 
(4,2,1) 

(4,2,1) 

4,2,1 
4,1,2 



(4,2,1) 

Then, as above, we can introduce four energy functions E^^, E^^, E'^^, E"^^ 
such that dE'^^ = u dv—x dy, dE^^ — u dv+y dx, dE^^ = —v du+y dx, dE"^^ — 
—V du — X dy (the superfixes correspond to the independent variables — thus 
for E^^ these are x and u , i.e., 1 and 3). 
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We will discuss these in more detail below where the rationale of our 
notation will be explained. In terms of the classical notation: 



dE = T dS — pdV (energy) 

dF = -SdT - pdV (free energy) 

dG = —SdT + V dp (Gibbs potential) 

dH ^ TdS + Vdp (enthalpy), 

i.e., E^^ = G, E'^^ = F, E^^ = H and E^^ = E. 

If we arrange the energy functions in lexicographic order , i.e., as E^^, 
E^^, E"^^, i?^^ and denote them by 5, 6, 7 and 8 in this order, then we can 
incorporate them into our system. For it follows from the definitions and 
simple substitutions that 

dE^^ = {y-vfi)dx- f2vdy 
dE^^ = (ufi + y)dx + uf2dy 
dE'^^ = -vfidx + {x - vf2)dy 
dE^'^ — ufidx + {uf2 — x)dy 

and so 

(5,1,2) = y-gh, (5,2,1) = -gf2 
(6,1,2) = y + fgi, (6,2,1) = fg2 
(7,1,2) = -gh, (7,2,1) = -x-gf2 
(8,1,2) = fg^, (8,2,1) = -x + fg2. 

One of the potentially irritating features of the thermodynamical iden- 
tities is that many are related by a simple swapping of the variables while 
this is accompanied by changes of sign which seem at first sight to be ran- 
dom. The simplest example is displayed by the four Maxwell relations. We 
can systemise such computations by introducing the symbol [a, 6; c, d\ for the 
Jacobi determinant of the mapping (c, d) i-^ (a, 6), i.e., 

[a, 6; c, (J\ — (a, c, d){h, d, c) — (a, d, c){b, c, d). 

The determinant then takes care of the sign. 

For example [3,4; 1,2] is the Jacobian 7—-^ — - and is therefore 1 (which 

d{x,y) 

Sill 1 f 
here denotes the number 1), [3, 2; 4,1] is ^, ' . and therefore = — ^ = 

d[v, x) g2 

(3) 1) 2) 

— ' ' . Note that there are 1,680 such Jacobians. However, lest the 

(4) 2, 1) 

reader despair, we then have the following simple rules for manipulating 
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these expressions which allow us to express them all in terms of our primitive 
quantities (/ and g together with their partials and, of course, x and y). 

[a, 6; c, (£\ = —[b, a; c, d] = —[a, b; d, c] 
[c, d] a, b] 



[a, b; c, d] 
(a, b, c) = [a, c; b, c] 
Further useful rules for computation are 

0,6; c,d\ = -r— 

[c, d; e, /] 

in particular, 

[a,b;e,b] 
[c, d; e, 6] 

and 

Using these rules, we can compute any of the 336 expressions of the form 
{i,j,k) (for i, j and k running from 1 to 8) by routine computations. For 
example, if we wish to compute (8, 3, 5) then we proceed as follows: Firstly, 

(8,3,5) ^[8,5;3,5] = 

Now [8, 5; 1, 2] = (8, 1, 2)(5, 2, l)-(8, 2, 1)(5, 1, 2) and [3, 5; 1, 2] = (3, 1, 2)(5, 2, 1) 
(3,2,1)(5,1,2). 
Hence, finally 

(oo.^ (8,1,2)(5,2,1)-(8,2,1)(5,1,2) 
^^""'^^ (3,1,2)(5,2,1)-(3,2,1)(5,1,2) 

and so can be expressed in terms of / and g, together with their partials. 
As a simple example, we can compute again the basic formulae 

(4,3,l) = f, (4,1,3) = (2,3,1) = ^, (2,1,3) = ^. 

J2 J2 J2 J2 

For example 

and the other three terms can be computed analogously. 
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5.2 Higher derivatives 



Some of the thermodynamical identities involve higher derivatives and we 
indicate briefly how to incorporate these into our scheme. We use the self- 
explanatory notation ((a, b, c)^d, e) for second derivatives. Thus ((3, 1, 2), 2, 1) 

is just /i2. Note that this notation allows for such derivatives as 

which is ((8, 1, 2), 3, 4). Once again, we can express all such derivatives (there 
are now 18,816 of them) in terms of x, /, g and their partials (now up to 
the second order) using the chain rule. For 

((a, 6, c),i,j) = ((a, b, c), 1, 2)(1, + ((o, b, c), 2, 1)(2, 

and (a,6, c), and (2,z, j) can be dealt with using the above tables. 



( - ( - ) ) 



5.3 Derived quantities and thermodynamical identi- 
ties 

The reason why there is a plethora of thermodynamical identities is simple. A 
large number of significant (and also insignificant) quantities can be expressed 
or defined as simple algebraic combinations of a very few (our primitive 
quantities .t, y. f, g and their partials). Hence there are bound to be many 
relationships between them. Our strategy to verify (or falsify) an identity is 
to use the above methods to express both sides in terms of these quantities 
and check whether they agree. 

Of course, there are myriads of such quantities and identities and we can 
only bring a sample. Thus we have 

the heat capacity at constant volume, and 



V 



dT J 



p 



the heat capacity at constant pressure. In our formalism, cy = /(4, 3, 2) and 
Cp — /(4, 3, 1), and so, from our tables, 

cy = /-p, Cp = f — . 

Hence for the important quantities 7 = — and Cp — cy we have 7 = 

^ cv im 

and Cp-cv^ f~rT- 
/1/2 
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Further examples are 



^v-(l^) =-(4,2,3), 

T 



dV 

the latent heat of volume increase, and 



the latent heat of pressure increase. 
Further definitions are: 



mv-{%] =(4,2,1) 



and 



V 

The coefficient of volume expansion at constant pressure is 

and the isothermal bulk modulus of elasticity is 

fdP\ 

Bt^-v[ — ] =-y(l,2,3). 



\dV 



T 



1 -1 

Then Kt — -=— — r is the isothermal compressibility. 

Bt Z/(1,2,3) 

We illustrate our method by verifying the simple identity: 

<^,-c.=r(g^g)^. 

Using the tables above, we can easily compute both sides in terms of our 

/ 

primitive expressions and get — — in each case. 

/1/2 

5.4 Computing (a, 6, c) and ((a, 6, c), c?, e) 
We can summarise these results in the following formulae: 

(a,l,2)(6,2,l)-(a,2,l)(b,l,2) 



[a, 6; c, (J\ 



(c,l,2)(d,2,l)-(c,2,l)(rf,l,2) 
21 



and so 

(a,l,2)(c,2,l)-(a,2,l)(c,l,2) 



(a, b, c) = [a, c; b, c] — 



(6,l,2)(c,2,l)-(6,2,l)(c,2,l) 



which allow us to systematically compute any of the derivatives of the form 
(a, b, c) in terms of our primitives x, y, f, g, fi, /2, gi and g2, using the above 
data basis for expressions of the form (a, 1, 2) and (a, 2, 1). 
For the second derivatives we substitute 

(a,l,2)(c,2,l)-(a,2,l)(c,l,2) 
(&,l,2)(c,2,l)-(6,2,l)(c,l,2) 

into the formula 

(0,l,2)(d,2,l)-(0,2,l)(d,l,2) 
'^^ (rf,l,2)(e,2,l)-(d,2,l)(e,l,2) 

to compute ((a, fe, c), d, e) in terms of our primitive terms (this time with the 
first and second derivatives of / and g). The advantage of these formulae is, 
of course, that one can write a simple programme to compute them. (It is 
always tacitly assumed in the above formulae that the appropriate conditions 
which allow a use of the inverse function theorem hold). 

It remains only to produce the corresponding data basis for second deriva- 
tives, i.e. to express all of the non-trivial quantities of the form ((o, b, c), d, e) 
with 6, c, d and e either 1 or 2 in terms of y and / and g and their partials. 
Of course, ((3, 1, 2), 1, 2), ((3, 2, 1), 1, 2), ((3, 2, 1), 1, 2), ((3, 2, 1), 2, 1) are just 
fill fi2 (twice) and f22- Similar identities hold for the partials of g. 

Further, 

((5, 1,2), 1,2)) = -gj,-gf,,; 
((5, 1,2), 2,1)) = l-g,f,-gf,,; 
((5, 2,1), 1,2)) = -gj,-gf,,., 

((5, 2,1), 2,1)) = -^2/2-^/22. 

Note that the two expressions for the mixed partial coincide, since fig2 — 
f2gi = 1. 



Similarly, 



((6, 1,2), 1,2)) = fm + fgu; 

((6, 1,2), 2,1)) = l + f,g, + fg,,- 

((6, 2,1), 1,2)) = fm + fgu; 

((6, 2,1), 2,1)) = g2f2 + fg22. 



22 



and 



((7, 1,2), 1,2)) = -gj,-gf,,- 
((7, 1,2), 2,1)) = -g,f,-gf,,- 
((7, 2,1), 1,2)) = -l-gj,-gf,,; 

((7, 2,1), 2,1)) = -^2/2-^/22. 

((8, 1,2), 1,2)) = g^h-fgn; 
((8, 1,2), 2,1)) = gif2 + fgi2; 

((8, 2,1), 1,2)) = -l + j\g, + fg,^; 

((8, 2,1), 2,1)) = -g2f2 + fg22. 

We emphasise that the numerical code for the various thermodynamical 

quantities is a mere construct to faciUtate their computation (ideally with 
the aid of suitable software) and that the final goal is to express them all in 
terms of the basic quantities (x, y, f, g and the partials of the latter). It 
is then a routine matter to translate these into the standard terminology of 
thermodynamics if so required. 



5.5 A notational survival kit 

In this treatment we have used three notations — the mathematically neutral 
symbols x, y, u and v etc. (to develop the mathematical theory which is 
independent of any reference to thermodynamics), the standard thermody- 
namical terminology p, V etc. (for readers interested in the thermodynamical 
interpretation) and finally the numerical code (to systematise the computa- 
tions of derived quantities and our approach to thermodynamical identities) . 
For the convenience of the reader we give a dictionary of the relationships 
between them: 

1 2 3 4 5 6 7 8 
p V T S G H F E 
X y u V E^^ E^^ E^^ E^\ 

Of course, p is pressure, V volume, T temperature, S entropy and G, H, F 
and E are free enthalpy, enthalpy, free energy and energy respectively. 



6 Existence and uniqueness of ^'-transversals 

In this section, we consider in more detail the restraints which are imposed on 
families of curves by the Samuelson area condition. We begin with a special 
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situation which we can compute directly and then show how to reduce the 
general case to it. We show that, given any family of level curves, there always 
exist families of intersecting level curves which satisfy the area condition, and 
the method of proof allows us to compute these curves explicitly in many 
interesting cases. Speaking loosely, this means that in thermodynamics to 
every family of isotherms (adiabats) there correspond families of possible 
adiabats (isotherms) (this subject is discussed more carefully below). Since 
the area condition is not very demanding, there are infinitely many collections 
of level curves which arc S-transvcrsal to a given family; but we show that if 
we are given two curves which are transversal to the latter (in the differential 
geometric sense) then they can be embedded in an essentially unique fashion 
into a system of S'-transvcrsal curves. This allows us, for example, to write 
down all families which arc S'-transvcrsal to the set of isotherms for the van 
der Waals gas. We include some examples which we found to be of interest 
in a later section. 

6.1 A special case 

We begin by investigating the questions of existence and uniqueness when 
one of the families consists of lines parallel to one of the axes, in this case, the 
X-axis. For this example, elementary computation shows that, as claimed, 
we can always describe all other possible families of level curves which are S'- 
transvcrsal to the first one. Moreover, in this case, it is also straightforward to 
show that knowledge of two curves determines the whole family. Remarkably, 
as we show in the next section, the general case can be reduced to this one, 
thus allowing a simple derivation of the basic theorems. If the first family 
of curves is calibrated, i.e., they are the level curves of a particular potential 
function u, then any single curve suffices to determine the second family. 

So let the f- foliation consist of the lines parallel to the x-axis (i.e., where 
v{x, y) = y), where in order to avoid topological difficulties we suppose that 
our potentials are defined on a product of intervals. (We are exchanging the 
roles of u and v here and will compute all w-functions which are transversal 
to this v). 

Then we know that if the foliation induced by the potential u is S'- 
transvcrsal, we can recalibrate u and v so that the Jacobian is identically 
one. If we assume that the u foliation has already been recalibrated and that 
the recalibration of v is v{x,y) — c{y), then a straightforward computation 
shows that J — c!{y)ux and so u must have the form: 

u{x,y) = a{y)x + b{y) 
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where b is an arbitrary smooth function of one variable and a is such that 

a(y)c'(y) = 1. 

We can state this formally as follows: 

Theorem 5 The function u is S\-transversal to the function v — c{y) if and 

only if u has the form u(x,y) = a(y)x + b(y) where aiy) = — -— . (Note that 

we are assuming here that c is a diffeomorphism between two intervals of the 
line). The function u is S -transversal to the function v = y if and only if u 
has the form 

u{x,y) = (f){a(y)x + b{y)) 

where (f) is a diffeomorphism between two intervals of the line, and a and b 
are any two smooth functions of one variable, for which a has no zeros. 

Since it will often be convenient to switch the roles of u and v, or x and 
y respectively, in the above, we document the corresponding formulae: 

u{x,y)^c{x), v{x,y) = a{x)y + b{x). 

We now consider the question of uniqueness in the above situation. By virtue 
of the general theory developed in the next section, this will suffice to cover 
the general case. Our starting point is the typical pair of S'l-transversal 
functions 

u{x,y) = a{y)x + b{y), v{x,y) = c{y) 

for arbitrary (generic) functions a and b (of one variable), with ac' = 1 (i.e., 

c is a primitive of — ) for the case where the i;-lines are the parallels to the 
a 

We now suppose that we have two transversals to the which we 

want to incorporate into a family of geometric adiabatics. We can suppose 
that the curves correspond to the values u — and u — 1. If c = 0, then we 
have 

and, for c = 1, 

^ _ 1 - bjy) 
«(?/) 

We note now that if the it-level curves are to be -S'l-transversal to the parallels 
to the X-axis, then they are transversal in the differential geometrical sense 
and so can be regarded as the graphs of functions (more precisely, 
function of y). Hence if we suppose that the "adiabatics" c = 0, c = 1 have 
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the form x = fo{y), and x = fi{y), then a simple computation shows that 
the general adiabatic has the form u = c, where 



6.2 A reduction 

We now show that the general case can be reduced to the previous special 
case. We begin with the remark that if we have any smooth non-vanishing 
function / of two variables, say on the product of two intervals, then we can 
always find a smooth vector field (m, v) which has / as its Jacobi function. 

Probably the easiest way to do this, as was pointed out to us by Michael 
Schmiickenschlager, is to use a field of the form 

u{x,y) = (t){x,y), v{x,y) = i/jiy) 

with a smooth function of two variables, ip one of one variable. (Such 
fields are called Knothe fields). The above form also has the advantage of 
leaving the level curves y = d invariant. The Jacobian of the above function 
is (j)x{x, y)il''{y) and we can, of course, easily choose the two free functions in 
such a way that this product gives /. For example, in the case where / is 
the constant function 1, then we can take for -0 any smooth function of one 
variable and then (j) is determined up to a function of y alone , i.e., has the 
form (p{x, y) = + x{y) where x is an arbitrary smooth function of one 

varia bl(0 

Using this result, we can prove the following: 

Theorem 6 Suppose that we have a foliation of part of the plane by the level 
curves of a suitable function u{x,y) (with non-vanishing gradient) which is 
defined on a domain (i.e., an open, connected subset) G in R^. Then we 
can linearise u locally by means of an area-preserving mapping. More pre- 
cisely, for each point (xo,yo) G we can find a neighbourhood G of the 
point in G and a function v{x,y) on G which is such that the mapping 
{x,y) I—)- {u{x,y),v{x,y)) is area-preserving and maps the lines u = c onto 
lines parallel to the y-axis. 



-'^When we include such a formula we are, of course, tacitly assuming that the operations 
carried out on the generic functions involved are legitimate. In this case this means 
explicitly that we are assuming that the derivative of "0 never vanishes , i.e., that "0 is a 
diffeomorphism. This type of situation will occur frequently in the following and since it 
would be tedious to state the explicit assumptions on the generic functions which arise, 
we will rely on the reader to fill in the details. 
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This is another result which is part of mathematical folklore. 

In order to prove it, we start with a foliation consisting of the level curves 
of a potential u and a point (xq, yo)- Since the gradient of u never vanishes, 
we can find a function v which is transversal to in a neighbourhood of this 
point and we introduce the new variables X = u{x, y), Y = v{x, y). By the 
inverse function theorem we can suppose that this can be solved to obtain x 
and y as smooth functions of X and F, say x = a{X, Y), y = b{X, Y). 

We now introduce further new variables X and Y of the form 

X = ^{X), Y = ij{X,Y) 

for suitable smooth functions and ip of one and two variables respectively. 
Then elementary calculations show that the Jacobian of X and Y with re- 
spect to the variables x and y (but expressed in terms of the variables X and 
Y) is 

<i>'{X)MX,Y)J{X,Y) 

where J{X, Y) is the Jacobian of {u, v) with respect to x and y, expressed 

as a function of X and Y via a and b , i.e., J{X,Y) = J{a{X,Y),b{X,Y)) 
where J{x,y) = Ux{x,y)vy{x,y) — Uy{x,y)vx{x,y). We can clearly arrange 
for this to be identically one by using the freedom in the choice of and i/j 
and this completes the proof. 

6.3 The general situation 

Using these results, we can now extend the existence and uniqueness results 
from the special case in which one of the foliations is parallel to the axis to 
any given it-foliation. 

Since the method is explicit we can also use it to find all possible S- 
transversal systems for several interesting special types of -u-curves. The 
method used is as follows: Suppose that we can find an area-preserving map- 
ping which maps the it-curves onto the lines parallel the the x-axis. Then we 
can transfer the above example to this situation. We remark that this is in 
a certain sense a rigorous justification for a ploy of Maxwell's, who argued 
from this special situation (for reasons of simplicity), assuming that his con- 
clusions then carried over to the general case (cf. the passage: "For the sake 
of the distinctness in the figure, I have supposed the substance to be partly 
in the liquid and partly in the gaseous state, so that the isothermal fines 
are horizontal, and easily distinguished from the adiabatic lines, which slope 
downward to the right. The investigation, however, is quite independent of 
any such restriction as to the nature of the working substance", [Ma], p. 
155). 
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We state this formally as a theorem: 



Theorem 7 Suppose that we are given a foliation of the plane which we take 
to be the level curves of a suitable function u. Then there exists (locally) a 
family of curves (the level curves of a potential v) which are S -transversal to 
the level curves ofu. Furthermore, given any two curves which are transversal 
to the level curves of u, then there exists a unique family of S -transversal 
curves which include the given two. 



6.4 Examples of the uniqueness and existence results 

We bring some explicit computations in connection with the question of the 
existence and uniqueness of .S-transversals to some simple cases. 



The ideal gas: We begin with the case of the adiabatics for the ideal gas. 

In this case we use the new variables X = xy, Y = log (xy'^) to 

7-1 

reduce to the simple case of transversals to the parallels to the coordinate 
axes. 

Then we have that two functions u and v where u is a recalibration of 
xy'^ are S'l-transversal if and only if they have the form 

u{x,y) = a ( — ^—-{Inx + jlny)] xy + b ( ( — ^—{Inx + jlny) 
V7-1 / VVT-l 

v{x,y) = c ( ^—{Inx + ^\ny) 
\7 - 1 

where c is a primitive of -. 

a 

Similarly, two functions u and v where w is a recalibration of xy are 
S'l-transversal if and only if they have the form 

u{x, y) — c{xy), v{x, y) = a{xy) i (In a; + 7 In y) ] + h{xy) 

V7- 1 / 

where c is again a primitive of — . 

a 

From these formulae it is easy to give the general form of functions u 
which are S-transversal to the adiabatics of the ideal gas resp. functions v 
which are S-transversal to its isotherms. 

At this point we bring a concrete example related to the ideal gas which 
was constructed to answer a question of Samuelson. Suppose that we are 
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given v{x, y) — xy and require an ^'-transversal function u which interpolates 
between the curves xy"^ = 1 and xy^ = 10 (i.e., two adiabatics corresponding 
to distinct cases of the ideal gas) . Then a simple computation shows that 

21n(V) 

is such that the contour w = is the first curve, while u — lis the second one. 
Interestingly, this then forces u to contain adiabatics of all the intermediary 
exponents, as the reader can easily verify. (This is an example where the 
representations are only valid locally, since any two curves of the form xy^ = c 
and xy^ — d will cross) . 

Analogous considerations lead to the following result: 

Theorem 8 Let v be S -transversal to u = xy. Then if two level curves of v 
have the form xy'^ constant for a fixed j, v is a recalibration of xy'^. 

The van der Waals gas: We now turn to the van der Waals equation. In 
order to simplify the notation, we use the following solutions of the Maxwell 
relationships: 



V — — 

7 



In this case we use the new variables 



X^(x + ^^]{y-l) 



and 

1 /, / 1 



y=^(^ln(^x + -j+7ln(y-l) 

to reduce to the simple case. 

Then we have that two functions u and v where v is a, recalibration of Y 
are S'l-transversal if and only if they have the form 

uix,y) = aiY)X + b{Y), v{x,y) = c{Y) 

where c is a primitive of -. 

a 
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Similarly, two functions u and v where w is a recalibration of xy are 
5'i -transversal if and only if they have the form 

u{x,y) = c{X), v{x,y) = a{X)Y + b{X) 

where c is again a primitive of -. 

a 

Prom these formulae it is again easy to give the general form of functions 
u which are S'-transversal to the adiabatics of the van der Waals gas, or 
functions v which are S'-transversal to its isotherms. 

Similar methods can be applied to the Feynman gas which are now dis- 
cussed. 

7 Five basic models 

We conclude by collecting some explicit computations for various gas models: 
beginning with the ideal gas. 

7.1 The ideal gas 

Here the recalibration is u = xy, v = (Inxy'^). In this case we can 

7 - 1 

explicitly compute the relationships, which are obtained by permuting the 
variables to get: 

u 1 , , 

X — —, V — niu-\-my] 

y 7-1 

y = —, v = mw — mx; 

X 7 — 1 

7— 1 —1 7— 1 7— 1 

y = e~'"x~ , y = e~'" x~ . 

The reader can check that the four Maxwell relations are indeed valid and 
thus compute the corresponding energy fuctions. We shall shortly describe 
a simpler and more systematic way of doing this. 

7.2 A generalisation of the ideal gas 

We now consider a simple generalisation of the ideal gas. This will presum- 
ably not describe any real gas (but see the remarks on Nernst's law below): 
we include it since the results are particularly transparent. The starting 
point is the equation 

u = x'^y^, V = x^y'^. 
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with a 7^ 1, 6 7^ 1 and ah — cd^ 0. The canonical recahbration is 

U ^ J_ ( ^ \ ^a(d-c-J+\) b{d-c-J+l) 

^\d-c- J+lJ ' ^ 

V - —( ^c(a-b=J+l) d{a-b-J+l) 

~ ^[a-b-j+ij'^ y 

where J = ad — be, and so, for the special case J = I, (which we can always 
achieve by means of a simple recahbration) , 



}i rj,a(d—c) yb{d—c) 



d — C 
1 

a — b 



^c{a-b) yd{a-b) 



We have included this example since the results have a pleasing simplicity 
and symmetry (see, in particular, the further computations below). The case 
of an ideal gas can be obtained by setting a = c = 1 and letting b tend to 
1 but there are some subtleties involved, as the presence of the logarithmic 
term in the recahbration of the ideal gas would suggest. 

7.3 The van der Waals gas 

The natural cahbrations are 

In this case the computations for computing E"^^ and E'^'^ can be carried 
out by hand and we get: 

u a 1 / ,N 1 1 

X r + ^, V ^h\[y -b)^ -Inii, 

y-b 7-1 

and 

Again this can be used to compute the two energies, but see below for a more 
systematic treatment which gives all such functions. 



7.4 The Feynman gas 

Here u = xy and v = xy"'^^'^^ for a function 7 of one variable. This example 
was introduced by Feynman (see [Fe]) to cope with the fact that, in a real 
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gas, the adiabatic index depends on temperature. This is a case where one 
genuinely requires the equation in [Col] to verify that it is a Samuelson 
configuration. This turns out to be the case and the recahbrations (which 
are not computed by Feynman) are 

U = (pixy) V = ln(xy^(^^)) 

where is a primitive of . This is the first example which we have 

7-1 

met where we genuinely have to recalibrate temperature (i.e., Boyle's law 
holds only in the weak form that pV is constant for constant temperature). 
For a discussion of the relevance of such recahbrations, see Chang [Ch] . The 
recahbrations introduced provide an at least qualitative explanation of the 
diagram on p. 78 of this reference, which displays comparative data of Le 
Due on spirit thermometers. 

Here we can solve for the cases where u and x, or u and v, are the 
independent variables. (We would like to thank P.F.X. Miiller who pointed 
out this passage in Feynman's text to us). 

7.5 A synthesis 

We can include all of the above (except the second example) in the form: 

u^(^x+^^{y-b) v^(^x + ^^{v- 6)^('^(^'^)) . 
with recahbrations 

u = (t>((x + ^]{y-b)] v = ln( (x + -^]{v- 6)^(<^"'("(^'^)) 



yV / vv 



Once again, is a primitive of . This is another case where it seems 

7-1 

hopeless to check that this represents a Samuelson configuration without the 
theory and computational methods developed here. 



7.6 A gallimaufry of formulae 

For completeness, we now bring a list of the expressions {i, 1, 2) and {i, 2, 1) 
for the substances introduced above. We emphasise again that we include 
these results since they allow us to compute the various energy functions 
and further thermodynamical quantities without explicitly calculating the 
various permutations of the variables implicit in the definitions. We have 
found no indication in the literature that this is possible. Again we start 
with the ideal gas: 
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The ideal gas: For the ideal gas 

(5,1,2 

(5,2,1 
(6,1,2 
(6,2,1 

(7,1,2 

(7,2,1 
(8,1,2 
(8,2,1 



ylog (xy^) 

y — ; 

7-1 
xlog {xy'^) 

7-1' 

7a; 

(ylog {xy^)) 

(7-1) ' 

{x\og{xy^)) 



-X 



(7-1) ' 



y 



(7-1)' 

X 



The generalisation of the ideal gas: Here 



(5,1,2)- 

(6,1,2) 
(7,1,2) 
(8,1,2) 



hy 
b — a 

dy 
d — c 

ay 

b — a 

cx 
d — c 



(5,2,1) = 

(6,2,1) 
(7,2,1) 
(8,2,1) 



bx 
b — a 

dx 
d — c 

ax 
b — a 

cy 
d — c 



12 



bxy 
b — a 



E 
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dxy 
d — c 



^23 



Hence the energy functions are given by E 

^ and E-^ = 
b — a d — c 

Note the pleasing symmetry of these results. We have found them useful 

as a litmus test for the validity of thermodynamical identities. 

We remark here that although this model may not correspond to any real 

gas, it does have the advantage that it satifies Nernst's law — the third law 

of thermodynamics in the precise form given in [La] — i.e., that if we express 

the entropy as a function of p and T or of F and T, then we get the form 

S — Pi(p)T" or = P2{p)T"^ for suitable positive indices n and m and 
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functions Pi and P2 of pressure. We know of no other explicit model for a 
real gas which has this property. 

A non-example Continuing on the theme of Ncrnst's law, we note that we 
have examined the Feynman model in this respect (for natural choices of 7) 
and found that it again failed to reproduce this phenomenon — the problem 
lies in the logarithm term in the recalibration of entropy. In view of the 
above remark, it was then tempting to combine the Feynman model and the 
above generalisation of the ideal gas, i.e. to consider the case 

Unfortunately, these functions do not normally satisfy the S'-condition. De- 
spite this disappointment, this computation at least shows the usefulness of 
the P.D.E. characterisation of the latter condition, in particular that it can 
be used to eliminate possible models which cannot be recalibrated to satisfy 
the Maxwell relations. 
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The van der Waals gas: 

-b + y) 



5,1,2) = 



5,2,1) = 



7,2,1) 
8,1,2) 



(7-1) ' 



y3 y2 J yyy2 



7-1 

{{-b + y)log({x+-^){-b + yr)]) 
6,1,2) = y ^ y ^ ■ 



(7-1) 



yd y2 J yyy2 



6.2.1) = - ^ 

7-1 

7.1.2) = y- 



7-1 



(.-^)--(7(^+^)(.-^r--g^) 

7-1 

?/ - ^ 
7-1' 

(y-6)^-^ (7 (^+^) (2/-^) 



'j.-l _ 2a{y-b) 
V 



1 

3 



8,2,1) = '--X. 

7-1 
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The Feynman gas: 



(5,1,2 
(5,2,1 

(6,1,2 

(6,2,1 

(7,1,2 
(7,2,1 

(8,2,1 
(8,2,1 



y-ylog {xy^^'^y'^) (^'(xy); 
-a; log (xy^^^^^) (j)'{xy); 

(f){xy) (y^^'^y^ +xlog{y)-f'{xy)y^^'=y'^+^) y-^^'^y^ 



x 



+ y; 



= (l){xy) 



j{xy) 



+ x\og{y)'y'{xy) 



-I/log {xy^^'^y^)(^'{xyy, 

— log (x?/'^^^^^) (f)'{xy)x — x; 

y-ii.^y)(j)^xy) + x\og{y)i{xy)y^^''y^+^) _ 



X 



(p{xy) 



l{xy) 



V 



+ x\og{y)i{xy) - x 



In reading Feynman's treatment, one gains the impression that he is tacitly 
assuming that the formulae for his model are obtained simply by plugging 
a variable 7 into those for the ideal gas. The presence of terms involving 
the derivative of 7 in the above show that this is not the case (for example, 
in the formulae for the important quantities Cp, cy and their difference and 
quotient). 

It is an easy task to compute the above quantities for the combined Feyn- 
man and van dcr Waals gas (using Mathematica) , but the results are too 
elaborate to be included here. 



8 Final remarks 

The mathematics of thermodynamics have never ceased to fascinate mathe- 
maticians, who generally experience a sense of unease at the standard repre- 
sentations, in particular of the laws of thermodynamics as an axiom system 
(see, for example, [Se] for a critical evaluation). There have been many at- 
tempts to put them on a solid basis. We mention, in particular, Caratheodory 
[Cal] and [Ca2], Lieb and Yngvason [Li] and Truesdell [Tr]. We have, of 
course, been influenced by these treatments and, inevitably, there are cer- 
tain common points. However, we believe that our approach is sufficiently 
original to justify its presentation. Thus in [Li] the ordering "adiabatically 
accessible from" is centre stage but apart from that the method is completely 
different. We know of two systematic approaches to thermodynamical identi- 
ties (Bridgman [Br] and Jayne [Ja]) and they have influenced our treatment. 
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Thus the idea of using Jacobians to derive identities can be found in the 
latteil. 

In conclusion, we would like t o express our gratitude to Iain Fraser and 
Elena Kartashova, who read and commented on an earlier version of our 
manuscript. 
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